STM conductance of Kondo impurities on open and structured surfaces 
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We study the scanning tunneling microscopy response for magnetic atoms on open and structured 
surfaces using Wilson's renormalization group. We observe Fano resonances associated with Kondo 
resonances and interference effects. For a magnetic atom in a quantum corral coupled to the confined 
surface states, and experimentally relevant parameters, we observe a large confinement induced effect 
not present in the experiments. These results suggest that the Kondo screening is dominated by 
the bulk electrons rather than the surface ones. 

PACS numbers: 72.15.Qm, 72.10.Fk 
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I. INTRODUCTION 

Scanning Tunneling MicroscopyB (STM) has proven to 
be one of the most useful tools to measure and charac- 
terize low energy electron spectroscopy, at least on open 
surfaces where the microscope tip has a direct access to 
the physical region undeij-, study. In a series of recent 
works, Kondo impuritiegd'S on the sur|£afje of noble met- 
als have been studied using the STM.QB'la In the Kondo 
problem there is a characteristic temperature Tx that 
separates the low temperature from the high tempera- 
ture regimes. In the low temperature regime {T < Tk), 
the spectral density of a Kondo impurity develops a res- 
onance at the Fermi energy, known as the Kondo reso- 
nance, and consequently locally changes the low energy 
spectroscopic properties of the system. When the micro- 
scope tip is placed on top of the impurity, the charac- 
teristic conductance vs voltage V shows a typical Fano- 
likeQia structure that is identified as the fingerprint of the 
Kondo effect. In this way, the Kondo behavior of several 
transition metal impurities on different noble metal sur- 
faces has been identified. In particular, the temperature 
dependence of the Kondo effect has been measured and, 
by studying Co impurities on different Cu surfaces, the 
scaling of the Kondo temperature with the host electron 
density has been shown Jj 

STM has also been used to create nanostructures on 
noble metal surfaces by manipulation of single atoms. 
This allows to positionate impurities at special points of 
the surface or to engineer the environment of the impu- 
rity. In recent experiments the Kondo effect for s^rGo 
impurity inside a quantum corral has been observed.E3 

The Fano like structure of the conductance observed 
when the STM tip is on top of a Kondo impurity, can be 
explained with a simple phenomenological model which 
includes an impurity resonant state at the Fermi energy. 
This resonance, that represents the Kondo resonance, has 
a width ksTK- In this model, the differential conduc- 
tance though the tip G {V) = dl/dV is given by: 



G{V) = go+a 



-1 + 2xq 
1 + 



(1) 



Here go is a background conductance, a is a constant. 



X = {eV — ef^I^bTk and q is the Fano parameter that 
depends on the electronic structure and on the tunneling 
matrix elements. 

For a more quantitative estimation, a many body cal- 
culation including the electron - electron interactions 
leading to the Kondo resonance is needed. Moreover, if 
the impurity is placed on a surface that is structured at a 
nanoscopic scale, the local density of states may strongly 
depend on energy, an—effect that changes the shape of 
the Kondo resonancetll and consequently of the voltage 
dependence of the STM conductance. A nanostructure 
with a characteristic length L has a characteristic en- 
ergy scale A ~ hvp/L where is the Fermi velocity. 
In systems with A ~ ksTx we expect strong effects. 
Quantum corrals on the Cu (111) surface can be engi- 
neered with the desired energy scale A, however it is 
not clear how much the surface states and bulk states 
contribute to the Kondo effect. Despite of some noto- 
rious effects like the observation, in elliptical corrals, of 
the so-called mirages,E3 some experiments indicate that 
the surface states play almost no role in the development 
of the Kondo effect.Q If so, the Kondo screening would 
not be affected by the corral that confines only surface 
states. A detailed analysis of the STM conductance can 
be used to elucidate this point, however to do that the 
surface nanostructure and the Kondo correlations have 
to be treated on the same footing. 

In what follows we present results for the tip conduc- 
tance on top of a Kondo impurity. Using the numerical 
renormalization group, we evaluate the finite tempera- 
ture conductance for an impurity on open and clean sur- 
faces and on structured ones. 



II. THE MODEL AND RESULTS 

In this section we present the model and briefly discuss 
how the STM conductance is calculated using Wilson's 
renormalization group. We then present the results for 
different situations. Our starting point is an extended 
Anderson model for magnetic impurities described with 
the following IIamiltonian:Ll 
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Ham = y^^egd^dg + Ud\d^^d\di ^^^£kc\^Cka 

+ E Vfcd\d„c\^,c^,,, (2) 

where the operator d\ creates an electron with spin cr at 
the impurity orbital with energy and Coulomb repul- 
sion \J , c! creates an electron in an extended state with 
quantum numbers fc and a and energy . The last term 
is a Falicov-Kimbal like term which takes into account a 
Coulomb repulsion between the electrons at the impurity 
orbital and the ones at extended states. The quantum 
number fc includes the band index and crystal momen- 
tum. In the case of impurities on the Cu (111) surface 
there is a surface band, that is relevant for determining 
the STM conductance.y 

The STM conductance at low temperatures is given 
byB 



(3) 



where eV is the voltage drop from the tip to the sample, 
and pt is the tip density of states that is assumed to be 
constant around the Fermi energy. In the linear response 
regime, the quantity pisp + eV) is given by 



p(eF + eV) 



1 



(4) 



Here Q{iij) is the retarded Green function of the oper- 
ator tcipti^) + tdd^r, tc and td are matrix elements for 
tunneling of an electron from the tip to the conduction 
band states and impurity orbital respectively and ip^{R) 
is the field operator that creates an electron on the sur- 
face at coordinate R. In what follows we consider that 
the tip is on top of the impurity and take R = Ri to be 
the impurity coordinate. 

The Green function is calculated using the numerical 
renormalization group (NRG). For a general description 
of the method we refer to Refs. |l^ and here we give 
only some details that are specific to our calculation. The 
starting point of the NRG is to reduce the Hamiltonian 
(^ to that of a linear chain 

Ham = Y.^ddUG + Ud]d^d\di+Y,Woflfoa 

+ ^ E(/olrf- + 4/0.) + Udc E dU.fL'foa' 
(J aa' 

+ E '^n(/no-/"+lcr + /rl+lcr/"<^) (5) 
(T,n— 

The states described by the operators are Wilson's 
orbitals centered at the impurity with the state |0, cr) as- 
sociated with /q^ being a Wannier like orbital centered 



at the impurity coordinate. The parameters A„ gener- 
ate a logarithmic discretization of the conduction band, 
and Wo describes the impurity potential scattering. The 
Hamiltonian (|^) is diagonalized numerically using an it- 
erative controlled scheme. At each iteration step iV, a 
truncated Hamiltonian with N orbitals is used and the 
many body states within an energy scale LOtq are evalu- 
ated. The thermodynamic and dynamical properties at 
each energy scale can then be calculated recursively. 

As it is usually done to describe the Kondo problem, 
Hamiltonian (|^) includes a single channel that is due to 
all the conduction bands. In the case of the Cu(lll) 
surface there is a surface band and bulk states that con- 
tribute to the density of states at the Fermi energy. Wil- 
son's orbitals are assumed to be a linear combination of 
surface and bulk states, in particular the Wannier state 
0,cr) contains a surface and a bulk component deter- 
mined by the respective hybridization matrix elements 

To calculate the STM conductance when the tip is on 
top of the impurity, we take the field operator ipl.{Ri) 
to be /q^ and calculate the retarded Green function of 
the operator t'^fQ^ + tdd^^. The matrix element t'^ depends 
on how much the Wannier orbital |0, a) contains of the 
surface or bulk states. To be more specific on the nature 
of the Wannier orbital lO,^) and on the matrix element 
t'^, a detailed model of the electronic structure of the 
metal host and of the electron transfer from the tip to 
the surface is required. 

The retarded Green function of Eq. (Q) is calculated 
as m Ref. |l| and 



p{eF + eV) 



1 



[Giu + iO+) - g{uj - iO+)]^^e^+eV 



A, A' 



Z{T) 



5{eF + eV - {Ey - Ex)) (6) 
where Z{T) is the grand partition function 

Z(T) = ^e-^-/'=-^ (7) 

A 

and A/a.a' = {X\t[fQ^ + tdd^^\X'). At each energy scale 
the discrete spectra is smoothed by replacing the delta 
functions by a smooth distribution. In what follows we 
use a logarithmic Gaussian that has been proven to give 
the best results. 



A. Impurity on a clean surface 

Here we present the results for a Kondo impurity on 
an open surface. In an energy scale of several ksTx , the 
bare density of conduction states is constant around the 
Fermi energy. Then the calculation can be done with the 
usual band discretization first proposed by Wilson. The 
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FIG. 1: STM conductance vs voltage, parameters in units of 
half the bandwidth are Sd = -U/2 = -0.5, Uds ^ 0, V = 0.2 
and (a) Wo = and different values of td/t'^ as indicated in 
the box; (b) td ~ and Wo ^ 0. In both (a) and (b) a Fano 
lineshape is obtained. 
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FIG. 2: Linear conductance spectra for different values of 
the impurity level energy Sd = —0.5U, —0.3U, —0.2U, and 
— 0.05C/. As Ed approaches the Fermi level, the Fano resonance 
becomes wider. Other parameters are (7=1, td/t'^ = 0.1, 
V = 0.2, and Wo = Uds = 0. 



results obtained for the conductance are shown in Figs. 
0, 1 and |. 

The STM conductance depends on the following pa- 
rameters: fcfiTif that gives the low energy scale where 
the local spectroscopic properties are dominated by the 
impurity, the temperature fc^T and the Fano parameter 
q that combines the tunneling matrix elements and some 
aspects of the band structure. If the Kondo resonance is 
represented by a simple resonant state, the Fano param- 
eter is given by 



-t^/mE(0) 



(8) 



where I](a;)— ,= J2k^k/{i^ — £k) is an effective band 
propagator.cJ For a system with electron-hole symme- 
try the real part of at the Fermi energy is zero and 
the parameter q is proportional to the matrix element 
td that describes the direct tunneling into the localized 
state. In this case, q accounts for the interference be- 
tween the two tunneling channels. In general, real sys- 
tems have a non zero ReY,{0) and there is no need to 
assume a direct tunneling to the localized state to have 
an asymmetric {q ^ 0) Fano resonance. In the NRG 
approach, the use of a symmetric conduction band sim- 
plifies the calculation and most of the numerical work 
on the Kondo effect has been done with this symmetry. 
We have introduced the potential scattering term Wq to 
break the electron-hole symmetry of the band. The Fano 
resonance then depends on the conduction band param- 
eters and in agreement with Ref. ^ with a non zero Wq 
we obtain asymmetric line shapes even for td — 0. 

In Fig. |l| the conductance vs applied voltage is shown 
for Ed — —U/2 and Udc — 0, all parameters are given 
in units of half the bandwidth D = 1. In Fig. |^(a) 
the results for Wq — and a set of the ratios td/t'^ are 
shown, this ratio determines the structure of the conduc- 
tance. In Fig. p](b), by independently varying Wq for 



= we show that in a full NRG many body calcula- 
tion, the low energy STM conductance is still given by 
a single parameter q as given by equation (|^). As Wq 
is increased the resonance narrows indicating a lowering 
of the Kondo temperature. This is due to a decrease in 
the local density of states at the Fermi energy. The in- 
crease of Wq also induces an increase of the Fano param- 
eter. For simplicity, in what follows we take symmetric 
band (i?el](0) — 0) where the Fano parameter is given 
hy q — td/t'^irpV, where p is the density of states of the 
band. 

Our calculation, that up to now mainly reproduces 
known results, shows that the NRG gives reliable re- 
sults when used to evaluate the STM conductance. In 
what follows we present some results that include the 
behavior of the conductance when the energy of the lo- 
calized orbital approaches the Fermi energy, the effect of 
the Falicov-Kimbal term Udc and the temperature depen- 
dence of the conductance. 

We start by discussing the conductance dependence on 
the localized orbital energy Sd- As Sd is increased and ap- 
proaches the Fermi energy, first the Kondo temperature 
increases and then the impurity enters an intermediate 
valence regime. In the conductance shown in Fig. |^ these 
effects are clearly observed. A notorious effect for inter- 
mediate values of Ed is the occurrence of an asymmetry 
for finite voltages: away from the Fermi level the values 
of G{V) are smaller for y < than for > 0. This ef- 
fect is due to the localized orbital resonance that is close 
to the Fermi energy. As Ed increases even more (not 
shown) the resonance crosses the Fermi level, the impu- 
rity occupation is reduced and the conductance shows a 
broad structure with a width that is given by the reso- 
nance width r rather than by the Kondo temperature. 
As we now show, the effect of a non zero electron-electron 
interaction Udc is equivalent to a shift in the impurity 
level. The parameter Udc is irrelevant, that means that 
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FIG. 3: Linear conductance spectra for different values of 
Udc = 0.05 (thin solid line), 0.1 (dotted line), 0.15 (dashed 
Une), 0.2 (solid thick line) and ea = -(7/2. (a) ta/t'^ = 0.1 
(b) td/t'^ ~ 0.3. Other parameters as in Fig. 2. 



in the Kondo regime the zero temperature fixed point 
does not depend on this parameter and the low tempera- 
ture thermodynamic behavior is universal as for Udc — 0. 
However this interaction may change some details in the 
low energy spectral structure for energies of the order or 
larger than kgTK- We have analyzed the effect of Udc on 
the STM conductance and the results are summarized in 
Fig. ^. As Udc increases the Fano line shape changes and 
two different effects can be observed: i) the width of the 
conductance minimum increases indicating an increase in 
the Kondo temperature and ii) the voltage dependence of 
the conductance becomes more asymmetric. These two 
points are consistent with an effective shift the localized 
level toward the Fermi energy. The asymmetry in G{V) 
obtained when the resonant level approaches the Fermi 
level has not been discussed in detail in the literature, 
however experimental results for Co on Cuflll) and for 
Ti on Ag show this type of effect. In Ref. the raw ex- 
perimental has a marked asymmetriy in the conductance 
that was attributed to a broad resonance lying just above 
the Fermi energy. Our results may explain the asymme- 
try as due to the same localized orbital that generates 
the Kondo resonance. 

Now we present finite temperature results and analyze 
the temperature dependence of the conductance. To do 
this we calculate the finite temperature spectral densities 
following Refs. |ll and |lj and use the finite temperature 
extension of Eq. 



G(T/,T) 



47re^ 



Pt I de 



dfje - eV) 
de 



p{e), (9) 



where /(e) is the Fermi function. In Fig. ^(a) the con- 
ductance for a small q and Wo — Udc = is shown for 
different temperatures. As the temperature increases the 
width of the low voltage structure increases. In agree- 
ment with the experimental results of Ref. |^, our results 
show a width that increases according to the low temper- 
ature Fermi liquid theory. In the inset of Fig. H (a) the 
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FIG. 4: (a) Linear conductance spectra for td/t'^ — 0.1 and 
different temperatures: T/Tk = (solid line), 0.75 (dotted 
line), and 1.5 (dashed line). Inset: Fermi liquid theory fit 
to the width 7 of the Fano resonance as a function of tem- 
perature, (b) Linear conductance spectra in an asymmetric 
situation Ed = —0.2, [7 = 1 and td/t'^ = 0.35 for different tem- 
peratures as indicated (the spectra at different temperatures 
have been shifted vertically for clarity). 



numerical width is shown for different temperatures, the 
thick continuous line is a fit with the Fermi liquid expres- 
^ = V2fcBTK Vi^keT)^ + 2{kBTKY. In Fig. |(b) 
the temperature dependence of a system with Sd = —0.2 
is shown. There is remarkable resemblance between some 
of these curves and the experimental results of Ti on Ag. 



B. Impurity on a structured surface 

We now present results for the STM conductance when 
the Kondo impurity is on a structured surface. The main 
effect of a nanostructure on the surface is to produce 
an energy dependence of the local density of states on 
a scale A that may be comparable with ksTx- If this 
happens the Kondo effect may be very much affected as 
discussed in previous works. Possible realizations of this 
situation may be obtained by placing impurities in quan- 
tum corrals, on small islands or in any structure with 
a characteristic size of the order of the hundred A. The 
problem then becomes much more complicated mainly 
because of two reasons: on the one hand, the nanostruc- 
ture of the host and the Kondo correlations have to be 
treated on the same footing. On the other hand, a de- 
tailed model for the metallic surface and the effect of the 
nanostructure on the bulk and surface states is needed. 
In general, in surfaces like the Cu(lll) where there is 
a surface band, we expect structures like the quantum 
corrals or small islands to confine surface states keeping 
the bulk states nearly unaffected. If so, these structures 
would change the Kondo effect according to the partici- 
pation of the surface states on the magnetic screening of 
the impurity. This point is still controversial: while the 



5 



0.8 



^0.6 



0.4 



0.2 





' 1 




- 











-0.0015 





eV 



0.0015 -0.0015 





eV 



0.0015 



FIG. 5: Linear conductance spectra for an impurity in a con- 
fined system Tk ~ A and S/A ^ 0.35 (a) at-resonance (b) 
off-resonance. Other parameter are Ed ~ —U/2, V — 0.2, 
td/t'c — 0.1, and Wo = Udc = 0. For comparison tfie line- 
shape for an open surface is shown ( thin line). 
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FIG. 6; Linear conductance spectra for an impurity in a con- 
fined system with td/t'^ = 0.1, Tk ~ A/10 and different 
widths 5 of the resonant states (a) at-resonance without bar- 
rier (thin solid line) and 5/ A ~ 0.63, 0.42, 0.21 (thick solid 
line) (b) off-resonance without barrier (thin solid line) and 
(5/A ~ 0.63, 0.42 (thick sohd line). 



observation of quantum mirages seems to indicate that 
surface modes actively participate in the Kondo screen- 
ing, other experimental resuhs suggest that the Kondo 
effect is due to bulk states. The STM conductance could 
be used to clarify this point. 

Due to the lack of a detailed model suitable to be 
treated with the renormalization group approach, we re- 
sort to a simple and phenomenological model to study the 
effect of structured local densities of states. The model 
reduces the problem to a single effective conduction band. 
As before, we assume that Wilson's orbitals are built as 
linear combination of the surface or bulk states. The sur- 
face nanostructure affects mostly the surface component 
of these states. In order to describe a local density of 
states with a structure in the desired energy range, we 
resort to the method described in Refs. |l^ and |l6|. The 
calculation includes a potential barrier that confines the 
electrons. The local density of states so obtained in the 
absence of impurity has the form of resonant states. This 
type of resonances have been observed using the STM. in 
atomic islands or terraces on the Cw(lll) surfaces.E3 If 
surface states participate actively on the Kondo screen- 
ing, the modulation of the relevant local density of states 
would be large indicating that the active electrons are 
confined by the nanostructure (a large potential barrier) . 
Conversely, if the Kondo effect is mostly due to the bulk 
states, the local density of states modulation would be 
very week. In this case we use a small effective barrier. 
The model is meant to make only qualitative estimations, 
nevertheless as we show below some important conclu- 
sions can be drawn from our analysis. 

As stated above, a nanostructure introduces a new 
length scale and the corresponding energy scale A. Since 
the number of microscopic parameters of the model is 
large, we just distinguish to cases: large structures with 
A < fcsTW and small nanostructures with A > /c^TV- 
Another important aspect of the Kondo effect in nanos- 



tructures concerns the position of the Fermi level relative 
to the structure of the local density of states. In this 
respect we consider the two extreme cases: the Fermi 
energy lying at a resonance (at-resonance case) or be- 
tween two resonances (off- resonance case) . The quantum 
corrals on Cm(111) with characteristic sizes of the order 
of lOOA correspond to the small nanostructure category. 
The elliptical corrals of Ref. ^ have the Fermi energy at 
a resonance. 

Figure ^ illustrates the case of a relatively large nanos- 
tructure (A ~ ksTK) for the at-resonance and off- 
resonance cases. The calculation was done with a small q 
and an intermediate barrier. The low temperature STM 
conductance shows a Fano like structure modulated by 
a superstructure of peaks separated approximately by 
the characteristic energy A. For the at-resonance case 
the central valley slightly broadens while for the off- 
resonance case narrow structures at low voltages are ob- 
tained. As the barrier height decreases the amplitude 
of the oscillations in the host local density of states de- 
creases and the conductance continuously approaches the 
open surface result. 

Figure || shows the results for a small system (A > 
/csTjf ). In this case the superstructure is not visible since 
A is out of scale. For the at-resonance (Fig. ||(b)) case 
and large barriers, the Kondo temperature aincreases and 
the Fano structure is broadened. This is due to an in- 
crease of the local density of states at the Fermi energy 
ep. The Kondo temperature has an exponential depen- 
dence with the density of states and consequently the 
broadening is very important for intermediate or large 
barriers. For the off-resonance case (Fig. ^(b)), the 
density of states decreases and the structure narrows, 
for large barrier a narrow Fano structure with a large q 
mounted on a background with a pronounced minimum 
at Ei? is obtained. 
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III. SUMMARY AND DISCUSSION 

We have analyzed the STM conductance G{V) for 
Kondo impurities on clean and structured surfaces. 
Starting with an Anderson model we used the NRG to 
calculate the physical properties and voltage dependence 
of the conductance. For clean surfaces we showed that 
the numerical approach gives reliable results and repro- 
duces the Fano type STM response. In particular, we 
have studied in some detail how G{V) depends on differ- 
ent microscopic parameters and on temperature. As the 
impurity level is shifted away from the symmetric case 
Sd = —U /2, either by changing its energy or by means of 
a Falicov-Kimbal interaction, the Fano line is distorted. 
For a set of physical parameters, both the shape and the 
temperature dependence of the conductance are in good 
agreement with experimental observations. The temper- 
ature dependence is also in good agreement with the low 
temperature Fermi liquid theory. 

To discuss the STM conductance of impurities in struc- 
tured surfaces we used an effective one channel model. 



In surfaces like the Cu(lll) one, where there is a surface 
band, this effective channel represents a linear combina- 
tion of surface and bulk sates. Since a surface structure 
perturbs much more the surface modes than the bulk 
modes, the sensitivity of the Kondo effect to the struc- 
ture gives an idea of the surface mode participation in the 
Kondo screening. We have shown that small structures 
like the quantum corrals, that are quite efhcient in confin- 
ing surface electrons, should produce important changes 
in the Kondo temperature and on the STM response if 
the surface modes were active in the Kondo physics. The 
experiments show that the Fano line has essentially the 
same width when the impurity is on a clean surface or 
inside a cpH:al, only small changes in the parameter q are 
observed.Ea This fact is a clear indication that the surface 
modes can sense the Kondo effect but are not important 
in determining it. This conclusion is in agreement with 
estimation oLthe microscopic parameters for surface and 
bulk states O 
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